CeBe macroscopically revisited 
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Magnetic susceptibility anisotropy and high sensitivity magnetostriction measurements are used 
to investigate the symmetry of CeBe ordered states. The antiferromagnetic state is confirmed as 
tetragonal, but no deviation from the cubic symmetry is observed in the so-called antiferroquadrupo- 
lar phase, where only volume effects are detected. In this phase, the temperature dependence of 
the strain field-susceptibilities is typical of non-ordered quadrupoles. Moreover, while an anti- 



ferroquadruplar order should be cubic, this symmetry is incompatible with the 



j) ordering 



wave-vector. The antiferroquadrupolar description of CeBg phase II is clearly inconsistent and an 
alternative model, based on a unidimensional representation of the cube, has to be sought for. 



In rare-earth compounds, the degrees of freedom of 
the 4/ ions are collectively involved at low temperature 
to form ordered states. The most common kind of order 
is the magnetic one that exists in a number of varieties. 
The orbital order is a less common situation, where only 
the orbital degeneracy of the 4/ ion is involved, with- 
out the stabilization of magnetic moments. It requires 
the orbital degeneracy to survive the crystal- field effect, 
which means the 4/ ions sites are of high symmetry. As 
the highest symmetry is that of the octaedra, cubic sys- 
tems are the most favorable scene for the development 
of orbital related effects There, the lifting of the 
degeneracy is equivalent to the emergence of 4/ elec- 
tric quadrupole moments and the associated order are 
referred to as " quadrupolar" . By analogy with the fer- 
romagnetism, the order that doesn't change the crystal 
periodicity is called ferroquadrupolar. One can also con- 
sider a state where the 4/ quadrupoles are changing from 
site to site, thus defining an "antiferroquadrupolar" or- 
der (AFQ), analog to the antiferromagnetism. For more 
than two decades, the most cited example of such an an- 
tiferroquadrupolar order has been CeBg phase II. 
This system crystallizes within the CaBg-type structure 
and displays complex magnetic phase diagrams [1, 
CeBg orders antiferromagnetically at T^v = 2.4 K (phase 
III) and its most intriguing feature is the occurrence of a 
non- magnetic state (phase II) below Tq= 3.2 K, before 
entering the paramagnetic state (phase I). In phase II, a 
magnetic field induced {W\) antiferromagnetism is de- 
tected via neutron diffraction [4| . Latter X-ray scattering 



experiments [5|-|8[ have confirmed a non-magnetic 
wave-vector, even in zero field. 

Inelastic neutron scattering results are consistent with 
a crystal- field splitting of the Ce'^"'", J = 5/2, multiplet 
resulting in a quadruplet ground state well separated 
(530 K) from the F^ doublet. As this ground state has 
an orbital degeneracy, Effantin et al. [J] interpreted the 
puzzling properties of phase II as resulting from an AFQ 
state. 

Up to now, no solidly based model has been given for the 



microscopic organization of this phase. Some authors 
suggest a Fg AFQ ordering while, for others, it is of 
the F^ type. The recurrent scheme is a F^ AFQ state 



where a single quadrupolar component, Pxy, is prop- 



agated by the [555] wave- vector. This model doesn't 
preserve the cubic symmetry and, due to the coupling 
of the quadrupoles with the lattice, a Fj^ ferroquadrupo- 
lar component and a macroscopic distortion should si- 
multaneously develop. Facing such inconsistencies and a 
confusing mass of experimental data, we thought that a 
crucial clarification would come from an accurate inves- 
tigation of phase II symmetry. 

This can be achieved from macroscopic experiments on 
a single crystal. For a rare-earth based system, mag- 
netic measurements are first to consider. In cubic crys- 
tals, an accurate determination of the first-order mag- 
netic susceptibility can reveal a symmetry lowering : 
any deviation from the cubic symmetry will result in 
an anisotropic susceptibility. Due to the magnetoelastic 
coupling, the lattice necessarily refiects a symmetry low- 
ering. The detection of the associated strain can be very 
sensitive, capacitance dilatometer being able to detect 
relative changes as small as 10~^. The difficulty, when 
investigating a sample undergoing a symmetry lowering, 
is that it divides into domains, mixing their properties. 
This can be solved using a magnetic field of appropri- 
ate direction, in order to coerce the sample to a single 
domain state. If the anisotropy of the susceptibility has 
been previously determined, the optimal field direction is 
easily defined. All the reported experiments where made 
using single crystals in the shape of rectangular platelets 
(6x4x1 mm'^) previously used for neutron scattering ex- 
periments [12|. 

The investigation of the magnetic susceptibility is based 
on magnetization measurements along the three main cu- 
bic directions ([001], [110] and [111]) in fields up to 8 T. 
The susceptibilities, xooij Xuo and Xm, are extracted 
using the Arrott method 13] applied to the descend- 
ing field data : only the field selected domains should 
contribute to the magnetization. Then, considering the 
tetragonal and trigonal symmetry lowerings [l3| cases, 
the measured susceptibility can be related to the normal 
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FIG. 1. Magnetic susceptibilities of CeBe along the three high 
symmetry directions of the cube, as deduced from Arrott's 
plots in descending field. 



susceptibilities X|| and x±y respectively parallel and 
perpendicular to the preserved, fourfold or threefold, axis 
of the system. One can thus identify the symmetry both 
via hierarchical and quantitative relations between xooi, 
Xiio and xm- 

Fig. [T] shows the thermal dependence of these three sus- 
ceptibilities at temperatures below 4.5 K. The transition 
between phases I and II is difficult to discern, both phases 
displaying a paramagnetic decrease of the susceptibility 
with the temperature. Within the precision of the mea- 
surements, typically 2%, no difference between the three 
susceptibilities and, therefore, no deviation from the cu- 
bic symmetry can be detected. In the antiferromagnetic 
state, the three susceptibilities are well separated, with 
a dominant Xooii an intermediate Xiio and a minimal 
Xiii- This sequence is expected in case of a tetrago- 
nal symmetry lowering, with a maximum susceptibility 
along the fourfold axis (x|| > x^)- In quantitative terms, 
this translates into xiio = \xooi + fxiiii which is pre- 
cisely what is observed on Fig. [TJ This is consistent with 
the multiaxial antiferromagnetic models Q with mag- 
netic moments perpendicular to the fourfold axis. 
Beyond the susceptibility analysis, one can further track 
a symmetry lowering using a high sensitivity dilatome- 
ter. For this purpose, we used the magnetostriction 
setup of the "Institut Neel", which is based on a ca- 
pacitance cell that can be rotated in a 0-6 T horizontal 
magnetic field, in the 2-300 K temperature range. The 
sensitivity, in terms of relative change in length, is bet- 
ter than 10^^. The alignment of the sample inside the 
setup is defined by the the choice of the crystal horizontal 
plane, which contains [mnp], the probed crystal's direc- 
tion, and [ijk] the direction of the applied field. For a 
reference length Iq , a current length / , the relative change 
in length, and its geometrical conditions, can be noted as 
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FIG. 2. Relative changes in length along [001] (a and b) and 
[111] (c) while rotating a constant field in the (110) plane in 
the ranges of phase III (a), phase II (b and c) and phase I (b 
and c). The zero angle a correspond to a field aligned with 
the measured crystal direction. The vertical lines indicate 
high symmetry directions in the (110) plane. 



: ™"PA,jfc = Al/lo = {l- lo)/lo- 

In the following, the default reference length Iq is taken 
at T = 5 K and in zero field. All the striction phenoma 
are analyzed in terms of cubic normal strain modes [l5| . 
For the most common cases of a tetragonal or trigonal 
symmetry lowering, the strain ™"pA of a single domain 
depends on three quantities: s", the volume strain, e'^, 
the tetragonal strain, and the trigonal one e^. 
To detect a symmetry lowering, one can take advantage 
of the adjustable field direction : rotating the sample in 
a constant field will successively select different domains, 
inducing changes in the probed sample length. Fig. [2] 
shows such measurements while rotating the field in the 
(110) plane, where all high symmetry axes of the cube are 
represented. In Fig. [2] (a) and[2](b), the probed sample 
direction is [001], which is ideal for detecting a tetragonal 
mode. Fig. [5] (a) shows what happens in the antiferro- 
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magnetic phase, under a 0.5 T applied field : spectacular, 
square like, changes in the sample length are observed. 
This is expected in case of a tetragonal symmetry lower- 
ing with maximum susceptibility along the fourfold axis 
: one goes from a single domain with [001] fourfold axis, 
for low OL values, to a a = ±90° state where the two other 
domains are selected, with fourfold axes [010] and [100]. 
The amplitude of the jump is directly related to the 
spontaneous strain of the antiferromagnetic state. Fig. 
[2] (b) gives the angular dependence under the maximum 
field (6 T) in phases II (T = 3.2 K) and phase I (T = 
6 K). These curves show no anomaly but a continuous, 
sine-type, change in the [001] length, which is typical of 
the paramagnetic state. The twisted phase II curve is an 
effect of higher-order susceptibilities in conjunction with 
a slight sample's misorientation. Fig. [5] (c) shows the 
angular variations, in phases II and I, of the [111] length, 
which is adapted to detect a trigonal mode. Here again, 
the continuous, undulating curves give no evidence of do- 
mains. Higher-order effects are clearly present in phase 
II, with a resulting double bump around a = 90°. 
The temperature dependence of the strain modes has 
been derived from a sequence of "^"'^Xijk measurements 
under a constant (0.5 T) field (Fig. |3]). It is only within 
the antiferromagnetic range that a spontaneous, tetrago- 
nal, strain is detected. The corresponding small accident 
on the trigonal curve is the result of an imperfect sam- 
ple alignment : a small fraction of the tetragonal strain 
is detected. In absence of a trigonal mode, the volume 
effect £" (inset of Fig. ^ can be safely derived from the 
length of a threefold axis : both transitions are visible, 
with a more pronounced thermal expansion in phase II, 
while phase III shows a substantial negative effect. This 
agrees with previous measurements by Sera et al. [l6| . 
The magnetic and magnetostriction data are consistent 
with a cubic symmetry in phases II and I. In this cu- 
bic context, the magnetostriction setup can be used to 
probe the magnetoelastic susceptibilities of the system 
with respect to a field enforced tetragonal (field along 
a fourfold axis) or trigonal (field along a threefold axis) 
symmetry lowering. This is the parastriction analysis 
of the paramagnetic state magnetoelastic properties [17| . 
At the lowest order, the strains vary quadratically with 
the applied field : = x^H^^ where is a magne- 
toelastic susceptibility and ^ states for 7 or e. In the 
harmonic treatment of 4/ magnetoelasticity, is pro- 
portional to the quadrupolar field-susceptibility Xq of ths 
same symmetry, (T) can be deduced by linearizing the 
constant temperature strain as function of the squared 
applied field. Fig. |4] shows examples of such curves in 
phase II (T — 2.7 K) and at the field crossing from phase 
I into phase II (T = 5 K). It appears that the latter has 
a surprisingly large volume dependence on the applied 
field (Fig. 21 upper part). On the e'^ parastriction curve, 
the transition from phase I to II is almost unnoticeable, 
whereas for the e mode, the negative slope is amplified in 
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FIG. 3. Temperature variation of the 7 and e strains, as 
obtained from measurements under a 0.5 T magnetic field. 
The inset shows the volume a strain deduced from the length 
of [111], in zero field. 



phase II at T = 5 K (Fig. |4l lower part). No hysteresis 
is observed in these curves, which confirms the absence 
of domain effects. In both phases, the absolute value of 
the susceptibility decreases with the temperature. This 
is expected if the degrees of freedom responsible for the 
field induced strains are not ordered. Quadrupoles are 
considered as the dominant actors in the parastriction of 
cubic rare-earth compounds, which results in a quadratic 
Curie-like behaviour for x'^(T). To test this behaviour, 
1 / y/x^ is represented as a function of T in Fig. \5\ In both 
phases, the curves seem indeed to vary linearly with the 
temperature. In phase II, the data for the 7 and e modes 
are almost superimposed. As these measurements mix 
elastic and magnetic properties, such an identity for dif- 
ferent symmetry modes is fortuitous. 

No symmetry lowering is detected in phase II, which 
discards 



non cubic models as the one of Ref. [11). 
While a true AFQ state needs to be cubic, the paras- 
triction results give evidence against the ordering of the 
quadrupoles. This is consistent with the elastic constants 
measurements [13], probes of the quadrupoles stress- 
susceptibilities, which, at least for the trigonal mode, be- 
have as if quadrupoles where not ordered. Nevertheless, 
we can deepen the analysis by specifying the character- 
istics of a cubic AFQ order. As the transition at Tg 
is second order, the ordered state should involve a single 
quadrupolar representation, Fg" (7) or F^" (e) , for a single 
wave-vector star. To obtain a high symmetry structure, 
stable down to K (above 2.5 Tesla), all the rare-earth 
ions should be equivalent by symmetry. 
Considering the Fj^ representation, the emergence of a 
4/ quadrupole is equivalent to a local tetragonal symme- 
try lowering : in a polar representation, the 4/ electronic 
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FIG. 4. Examples of parastriction curves in phases I and II, 
at T — 2.7 and 5 K, measured for the 7 and e symmetry 
lowering modes. Upper part : squared field dependence of 
the volume mode. Lower part : the tetragonal and trigonal 
strains. 
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FIG. 5. Thermal dependencies of l/y^X^ and l/y'x^ in the 
range of phase I an II. The lines are just guides to the eye 
and the inset details the data for phase II. 



distribution is stretched (or squeezed) along one of the 
fourfold axes of the cube. There are then three sym- 
metry equivalent Vt, modes for the 4/ ion. As for the 
representation, the quadrupoles develop with a local 
trigonal symmetry lowering, which means that there are 
four possible 4/ modes, one for each threefold axis. The 



cubic symmetry is preserved if all the modes, three for 
or four for T^, are equally represented in the cubic 
AFQ cell. The smallest cell one can consider is built from 
the 8 sites at the vertexes of a cube, which is enough to 
describe structures with (^00), (^^0) and (5^5) wave- 
vectors. On this cell, it is clearly impossible to equally 
represent the three Tj]' modes : there exists no cubic Vt, 
AFQ structure for these wave-vectors. This is not the 
case for and cubic AFQ states are possible for the 
(iOO) and (iiO) wave- vectors. However, this is strictly 
impossible for (^^^), the only wave- vector that is de- 
tected in phase II. This star reduces to a single branch 
which defines only two 4/ modes, with opposite sign. An 
object that alternates in sign from site to site, while pre- 
serving the cubic symmetry, necessarily belongs to a uni- 
dimensional representation of the cube. Among the even 
Oh group representations, only FJ allows a change in 
sign. To find an electric multipole transforming accord- 
ing to FJ, one has to go up to 64 poles, which means 
a sixth-order operator in 3 components. For J = 5/2, 
these hexacontatetrapoles cancel : no 4/ electric multi- 
pole can describe the cubic (^^^) ordering of CeBg. 
To solve this puzzle, one has to consider alternative ob- 
jects, transforming according to one of the unidimen- 
sional representations of Oh- In a cage system, as sug- 
gested by Kasuya [3], one may consider the dynamic 
decentering of the 4/ ion, already present above Tq as a 
Jahn- Teller effect. Below Tq there might be a symmetry 
breaking in the spatial distribution of the Ce ion, that 
would alternate from site to site. In contrast with the 4/ 
electronic configuration, the related electric multipoles 
can be of odd parity : in addition to Fj, they may also 
transform according to the or representations. 
The authors are indebted to Dr L.P. Regnault, who 
provided us with the sample and critically read the 
manuscript. 
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